Formation of the Q ball in the thermal logarithmic potential and its properties 
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We investigate the Q-ball formation in the thermal logarithmic potential by means of the lattice 
simulation, and reconfirm qualitatively the relation between Q-ball charge and the amplitude of 
the Affleck-Dine field at the onset of its oscillation. We find time dependence of some properties 
of the Q ball, such as its size and the field value at its center. Since the thermal logarithmic 
potential decreases as the temperature falls down, the gravity-mediation potential will affect the 
properties of the Q ball. Even in the case when the gravity-mediation potential alone does not allow 
Q-ball solution, we find the transformation from the thick-wall type of the Q ball to the thin-wall 
type, contrary to the naive expectation that the Q balls will be destroyed immediately when the 
gravity-mediation potential becomes dominant at the center of the Q ball. 



I. INTRODUCTION 



Q-ball formation is ubiquitous in the Affleck-Dine 
mechanism for baryogenesis llHal- Soon after the Affleck- 
Dine field begins rotation in the potential which al- 
lows the Q-ball solution, the homogeneous field starts 
to fluctuate and transforms into lumps. The actual for- 
mation was investigated numerically on the lattices for 
the gauge- and gravity-mediated supersymmetry (SUSY) 
breaking scenarios [JQ. 

On the other hand, the properties of the Q ball are 
quite different among the different forms of the poten- 
tial. For example, for the flat potential, such as in the 
gauge-mediated SUSY breaking scenario, the Q-ball size, 
the field value at its center, and the field rotation speed 
depend on the charge Q of the Q ball nontrivially p} . 

The potential may be dominated by thermal effects af- 
ter inflation. In the Affleck-Dine scenario, the field am- 
plitude at the onset of the rotation is very large, and 
the two-loop thermal effects on the potential are crucial 
[1, Q . This potential is given by 



(1) 



for large field values 1 . We called it the thermal loga- 
rithmic potential Q • The energy density of the universe 
is dominated by the oscillation of the inflaton after in- 
flation, but there exists dilute plasma [T2[ which would 
build up the thermal logarithmic potential. We consid- 
ered in Ref. [B[ that the Q-ball formation in the ther- 
mal logarithmic potential would be more or less similar 
to the time-independent logarithmic potential because of 
the fast growth of the fluctuations of the field in spite of 
the time dependence of the temperature T. We thus bor- 
rowed the results from the time-independent logarithmic 
potential case. For example, the charge of the produced 



Q ball is given by @ 



(2) 



1 We do not consider a negative thermal logarithmic potential llOi 
llll . because Q-ball formation does not occur in that potential. 



where fi rj 6 x 1CU 4 , and the subscript "osc" denotes the 
values of the variables at the onset of the oscillation of 
the field. 

In this article we actually perform lattice simulations 
of the Q-ball formation in the thermal logarithmic poten- 
tial, and see the evolution of the field and the distribution 
of the produced Q balls. We also study of the properties 
of the formed Q balls, especially about the evolution of 
the field amplitude at the center of the Q ball. 

In addition to the thermal logarithmic term, there 
would also be a mass term due to the gravity-mediated 
SUSY breaking effects for both the gauge- and gravity- 
mediation scenarios. The mass term alone (including 
one-loop potential) does not allow the Q-ball solution in 
some cases. In such cases, one may naively consider those 
Q balls, created when the thermal logarithmic potential 
dominates, to disappear when the mass term begins to 
dominate the potential at the field value of the Q-ball 
center. As shown below, however, the Q-ball solution 
still exists later. Actually, the thick-wall type (gauge- 
mediation type) Q ball transforms to the thin-wall type. 

The structure of the article is as follows. In the next 
section, we show the results of lattice simulations for 
both time-independent and thermal logarithmic poten- 
tials. Some properties of the Q ball with thermal loga- 
rithmic potential are shown in Sec. Ill, while, in Sec. IV, 
we focus on the transition from the gauge-mediation type 
to the thin-wall type Q balls when the mass term in the 
potential gradually dominates over the thermal logarith- 
mic potential. We finally conclude in Sec.V. 



II. Q-BALL FORMATION 

We investigate the Q-ball formation by means of three- 
dimensional lattice simulations. Interpolating the ther- 
mal mass term at smaller field amplitudes and the ther- 
mal two-loop potential growing logarithmically at larger 
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field values, we take the following form of the potential: 

y T ($)=T 4 log(l + ff). (3) 

In addition, we restudy the time-independent case where 
T is replaced by the constant mass Mp for comparison, 
which is nothing but the potential of the gauge-mediated 
SUSY breaking effects. Here we consider the inflaton- 
oscillation dominated universe after inflation before re- 
heating, but a similar argument applies also to the ra- 
diation dominated era. The temperature decreases as 
the universe expands as T oc a -3 / 8 oc t^ 1 / 4 , where a(t) 
is the scale factor of the universe. Since we are inter- 
ested in the period after the Affleck-Dine field starts the 
oscillation (rotation) when H 2 ~ \V"\, we take initial 
conditions as 2 

pi(0) = ¥Jo(l + *i), ip' 1 (0) = S 2 , 
¥> 2 (0)=* 3 , ^(0) = V2(l + <5 4 ) ! 

*°> = shf«" (4) 

where all the variables are normalized by the temper- 
ature at the onset of the oscillation, T osc , such that 
(p = $/\/2T osc , h = H/T osc , £ = T osc x and r = T osc t. 
Here we decompose the field into real and imaginary 
parts as ip — tpi + itp%. <5's represent the fluctuations 
of O(10 -7 ). We mostly use 256 3 lattices, but in order to 
see any box size effects we also perform on 350 3 lattices in 
some cases, but find no crucial differences between them. 

Figure [1] shows the initial amplitude dependence of 
the largest charge of the Q ball produced for the time- 
independent (lower) and thermal (upper) logarithmic po- 
tentials. Here we average the charge over five realiza- 
tions of the initial fluctuations except for the smallest 
initial amplitude case (<po = 300). The lower line cor- 
responds to the relation obtained in Ref. Q = (3ip 4 
where (3 w 6 x 10~ 4 , thus we reconfirmed the previous 
results. 

On the other hand, in the thermal logarithmic poten- 
tial, we have a similar relation Q = P'tp^ with j3' as 
2 x 10~ 3 , although it might have a little tilt. Thus, 
qualitative features for the Q balls in the thermal log- 
arithmic potential can be captured by the case with the 
time-independent logarithmic potential with Mp being 
replaced by T osc . 

The difference between the time-independent and ther- 
mal logarithmic potentials can be qualitatively consid- 
ered as follows. Q-ball charge can be estimated as 
Q ~ ql\ at the formation time, where q is the charge 
density and £h ~ t is the horizon scale. From Fig. [5J 
the formation times are do ~ 4.6 and ay ~ 9.3, re- 
spectively, for the time-independent and thermal loga- 
rithmic potential. Thus, the ratio of the charges would 



2 It is implicitly assumed that the helical motion is dynamically 
achieved by so-called A terms. 
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FIG. 1: Charge of the Q ball depending on the amplitude of 
the Affleck-Dine field at the onset of the oscillation. The lower 
(upper) points and line correspond to the time-independent 
(thermal) logarithmic potential. 

be ~ {a T /ao)- 3 (t T /to) 3 ~ (a T /ao) 3/2 ~ 3. Notice that 
the later rise and the slower growth of the amplitude 
of the fluctuations are due to the shrinking instability 
band < k/a < 2T 2 /c/>o and the lowering growth rate 
T 2 /(-\/20o) due to the decreasing temperature T. 
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FIG. 2: Evolution of the fluctuations of the Affleck-Dine in the 
time-independent (green, the earlier rise) and thermal (red, 
the later rise) logarithmic potentials for ipo = 10 s . Notice 
that all the parameters including initial conditions are taken 
to be the same. 

The distribution of the Q balls in example cases are 
shown in Table U for the thermal logarithmic potential 
and in Table |H] for the time-independent logarithmic po- 
tential, where we identified Q balls with Q > 10 3 and 
Q > 10 5 , respectively, in these cases. As one can see, the 
charge is dominated by those Q balls with the charge of 
the largest magnitude. Therefore, it is fairly reasonable 
to estimate any relations among the Q-ball parameters 
by using the largest Q ball, as we derived the charge and 
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the initial field amplitude above. Since it is beyond the 
scope of the present paper to provide an analytical esti- 
mate of the distribution, we just provide a fitting formula 
of the form 



Ni(Q) = OitpoQ-ie V 77 ^) (t = 0,T), (5) 

where Q denotes the charge in terms of the order of mag- 
nitude. Here r\ w 0.3. Q max = Pi'-Po with (3o w 6 x 1CT 4 
and Pt = 2 x 10 -3 for the time-independent and ther- 
mal logarithmic potentials, respectively, ao ~ 1.6 and 
ax ~ 1.1- 

TABLE I: Distribution of Q balls in the thermal logarithmic 
potential for 950 = 10 3 . 

Charge Numbers Sum of the charge Fraction 



O(10 9 ) 


2 


3.77 x 10 9 


0.7358 


O(10 8 ) 


2 


6.33 x 10 8 


0.1235 


O(10 7 ) 


15 


4.39 x 10 8 


0.0857 


O(10 6 ) 


15 


5.82 x 10 7 


0.0114 


O(10 5 ) 


12 


7.21 x 10 6 


0.0014 


O(10 4 ) 


51 


1.05 x 10 6 


0.0002 


O(10 3 ) 


195 


8.94 x 10 6 


0.0002 



where 

Vl ~THog(l + ^) ~T\ (9) 

V 2 ~uo 2 <j> 2 . (10) 

When the energy takes the minimum value, the equipar- 
tition is achieved for gauge- mediation type Q balls Q. 
From -Egrad ~ Vi ~ V2, together with the charge con- 
servation §§§ , we obtain the evolution of the (comoving) 
Q-ball size R, the rotation speed of the field to, and the 
field amplitude at the center of the Q ball tj> c , respectively 
as 

R cx a -(4-7)/(^+D ; (n) 

U CX -^-3+7)/(U+l) ) (12 ) 

^ cx a ^-i-^)/(D+i) ) (13) 

where we define 7 by Vi ~ T 4 cx a 1 . Notice that 
7 = 3/2 and 4 for the inflaton-oscillation and radiation 
dominated universe, respectively. These properties are 
observed in the lattice simulations which we perform for 
the D = 3 case. 



TABLE II: Distribution of Q balls in the time independent 
logarithmic potential for ipo = 3x 10 3 . 



Charge 


Numbers 


Sum of the charge 


Fraction 


O(10 10 ) 


5 


1.42 x 10 11 


0.7807 


O(10 9 ) 


10 


2.33 x 10 10 


0.1280 


O(10 8 ) 


20 


5.92 x 10 9 


0.0326 


O(10 7 ) 


44 


1.56 x 10 9 


0.0086 


O(10 6 ) 


56 


2.27 x 10 s 


0.0012 


O(10 5 ) 


163 


6.46 x 10 7 


0.0004 



III. Q-BALL PROPERTIES 

Let us investigate the evolution of the D-dimensional 
Q ball formed in the thermal logarithmic potential, where 
for D = 1 and 2 it is wall- and stringlike objects, respec- 
tively [3j . The charge of the Q ball is given by 



Q ~ a A R q ~ const., 



(6) 



where q is the charge density and R is the Q-ball size. 
Charge conservation implies that Q is constant. If we 



write $(x, t) 
written as 



E = I d 3 x 



= / d 3 x [E, 



1 



- Vi + V 2 ] + ujQ, 



(7) 



IV. TRANSFORMATION OF Q-BALL TYPES 

In addition to the thermal logarithmic term in the 
potential, there is a mass term which stems from the 
gravity-mediated SUSY breaking effects. This potential 
can be written as 



K\og 



(14) 



where m 3 / 2 is the gravitino mass, and one- loop effects 
are included. K is either a positive or negative constant 
of 0(0.01 — 0.1), and M* is a normalization scale. This 
potential alone allows the Q-ball solution only if K < 
[HEll]. In the opposite case (K > 0), Q-ball formation 
is prohibited. One may thus be apt to consider that the 
Q ball created in the thermal logarithmic potential will 
disappear once the mass term with K > dominates over 
the thermal one at the field value of the Q-ball center, 
Vt{4>c) < V m (<pc)- This condition can be written as 



> 



/ eq 



(x)e laJt /v / 2, the energy of the Q ball is Since 



3 / 4 and 



oc 



T 2 
WI3/2' 

a -3/8 f or d 



(15) 



= 3 in the 

inflaton-oscillation domination, and <fi eq cx a~ 2 and <p c cx 
a -1 for any D in the radiation domination, it is true that 
V m (4> c ) will eventually overcome Vt(4>c)- 

However, the Q-ball solution does exist for V = Vr + 
V m with K > for T > m 3 / 2 - Therefore, Q balls are not 
destroyed, but the metamorphosis will take place in such 
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situations. It might be best to simulate on the lattices to 
verify this phenomenon, but it is very time-consuming to 
perform. Here, instead, we take another approach, and 
leave the lattice simulations for future work. 

We seek the Q-ball solution for V — Vr + V m at some 
time snapshots. In order to obtain the solution, we just 
have to solve the equation 



^V + 2# + / 2 , 
dr 2 r dr 



(16) 



[U 



with boundary conditions 4>(oo) — and </>'(0) 

Since we would like to compare the results also to the 
solution of the pure logarithmic potential case, we change 
the value of the mass 7713/2 in V m , leaving the logarithmic 
potential time independent. Thus we use the following 
potential with various values for to 3 / 2 : 



V = Mp log 1 



2Mp 



1 



& 2M 2 J 
(17) 

Notice that it is practically the same if one varies 
Mp while TO3/2 is fixed, since the results are derived 
and shown by variables normalized with respect to 
Mp. The profiles are shown in Fig. [3] for m 3 / 2 /Mp = 
(0, 0.5, 1, 2, 3, 4, 5, 10, 20) x 10~ 7 from the top to the bot- 
tom. Here we set K = 0.1 and M^/Mp = 10 6 . It mimics 
the time evolution that the mass term eventually dom- 
inates over the logarithmic potential. Since the charge 
should be conserved, the angular velocity uj increases: 
oj/Mp = (2.0, 2.1, 2.4, 3.2, 4.1, 5.2, 6.2, 11, 21) x KT 6 from 
the top to the bottom. One can see that the thick-wall 
type (gauge-mediation type) of Q ball transforms into the 
thin-wall-like type as time goes on. Notice that 4> c > <fi eq 
takes place for m/Mp > 10 -7 . 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

M F r/10 6 



In the gravity-mediation, Q-ball solutions may exist 
until T ~ rn 3 / 2 for large enough uj. 3 After that, Q 
balls will disappear quickly, but the Affleck-Dine field 
can no longer be regarded as homogeneous because 
the field is localized near the place where destroyed Q 
balls had existed; they were very much separated from 
each other. Typical separation length is estimated as 
^_f/,formation(adcstruction/aformation)- On the other hand, 
in the gauge-mediation, Q balls will remain intact for 
77i0 > 777,3/2, where — \JV" (0) is the curvature of 
the gauge- mediation potential at = 0, since the Q-ball 
solution will exist irrespective of the temperature. 



V. CONCLUSIONS 



We have investigated the Q-ball formation in the ther- 
mal logarithmic potential by means of three dimensional 
lattice simulations. First of all, Q balls are actually 
formed. This is because the growth of the field fluctu- 
ations is fast enough to create Q balls, in spite of the 
shrinking instability band due to the decreasing temper- 
ature, and so on. We have found that the charge of the Q 
ball in the thermal logarithmic potential has almost the 
same dependence on the initial amplitude of the Affleck- 
Dine field, 



Q = P 



in it 



(18) 



with (3 s» 2 x 10 3 , which is a factor of 3 larger than that 
in the time-independent logarithmic potential case. 

We have also estimated the evolutions of parameters, 
such as the Q-ball size R, the field rotation velocity w, 
and the field value at the center of the Q ball <f> c . Since 
the thermal logarithmic potential decreases as the tem- 
perature drops, the mass term would dominate over the 
thermal logarithmic one at the field value <fi — 4> c . Even 
if the mass term V m with positive K alone does not al- 
low any Q-ball solution, the total potential V = Vr + V m 
does allow the solution. In fact, we have found that the 
thick-wall type Q ball will eventually transform into the 
thin-wall type. Finally, the Q balls will disappear when 
T ~ 777,3/2, with an inhomogeneous Affleck-Dine field be- 
ing left afterward in the gravity-mediation scenario. 



FIG. 3: Profiles of the Q balls with Q » 7.1 x 10 26 for 
m/M F = (0,0.5,1,2,3,4,5,10,20) x 10" 7 from the top to 
the bottom. The Q-ball profile "evolves" from the top to the 
bottom. 
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3 The fate of the Q ball will depend on the cosmological situation. 
For thorough analysis, see Ref. 
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